Dirac fermions at the H point of graphite: Magneto-transmission studies 
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We report on far infrared (FIR) magneto-transmission measurements on a thin graphite sample 
prepared by exfoliation of highly oriented pyrolytic graphite (HOPG) . In magnetic field, absorption 
lines exhibiting a blue-shift proportional to \/B are observed. This is a fingerprint for massless Dirac 
holes at the H point in bulk graphite. The Fermi velocity is found to be c = (1.02 ±0.02) x 10^ m/s 
and the pseudogap | A| at the H point is estimated to be below 10 meV. Although the holes behave 
to a first approximation as a strictly 2D gas of Dirac fermions, the full 3D band structure has to be 
taken into account to explain all the observed spectral features. 

PACS numbers: 71.70.Di, 76.40.+b, 78.30.-j, 81.05. Uw 



The fabrication of graphene, a 2D lattice of carbon 
atoms with a honeycomb symmetry, and the subsequent 
discovery of Dirac fermions P, 0, H, 01 > \ea,A to re- 
newed interest in the physical properties of bulk graphite. 
Compared to graphene, graphite represents a system of a 
higher complexity which is still not fully understood de- 
spite fifty years of intensive research. The currently ac- 
cepted tight binding (TB) model, formulated for graphite 
by Slonzewski, Weiss and McClure (SWM) 0,11], implies 
seven TB parameters 70, ... 75, A and predicts the pres- 
ence of particles with parabolic and linear dispersions at 
the K and H points of the Brillouin zone respectively. 
Whereas the appearance of massive electrons at the K 
point has been reported in numerous experiments, see 
e.g. 0, B HIj there is httle direct evidence for Dirac 
fermions (holes) at the H point. 

Evidence for -s/S-dependent features in magneto- 
reflectance spectra, typical for Dirac particles, has been 
reported by Toy et al. [l^. Recently, quantum oscilla- 
tions in the magnetoresistance of bulk graphite 
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indicated the presence of strictly 2D gases of both Dirac 
holes and massive electrons. Nevertheless, this inter- 
pretation remains controversial l^, 14 1 and raises some 
doubts concerning the accuracy of the SWM model. Fur- 
ther evidence of Dirac fermions close to the H point was 
obtained from angle resolved photoemission spectroscopy 
(ARPES) [ll,[l3- However, the ARPES measurements 
gave differing Fermi velocities, information on the band 
structure only below the Fermi level, and are charac- 
terized by a rather low accuracy. Dirac holes have re- 
cently been observed in scanning tunneling spectroscopy 
(STS) [13, indicating again the pronounced 2D behavior 
of particles in HOPG. However, these investigations are 
principally limited to the surface of the sample and it is 
not understood why these effects have not been observed 
in equivalent experiments performed earlier 
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In this Letter, we use FIR magneto-transmission exper- 
iments to probe the nature of the holes at the H point 
of HOPG. We address the currently controversial issue 
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FIG. 1: A Raman spectrum of the studied sample showing the 
Raman bands G and 2D typical of bulk graphite [l^. The 
bright laser spot in the middle of the microscope image in the 
inset shows the place, where this spectrum was collected. 



of whether the holes in HOPG can be described using 
a simplified model assuming a strictly 2D gas of Dirac 
fermions as proposed in Refs. [ll|, [l^, [l3|, or whether a 
full 3D band structure needs to be employed. We show 
that the appealing 2D model is in the simplest approach 
applicable, but simultaneously, we demonstrate the clear 
limits of this approximation. In this work, we focus on 
transitions, the energies of which scale as ^fB. However, 
in a different spectral region (low energies and high mag- 
netic fields) we clearly observe features linear in B, which 
in accordance to previous reports arise from electronic 
transitions in the vicinity of the K point. 

A thin sample for the transmission measurements was 
prepared by exfoliation of HOPG. A confocal microscope 
image of part of the sample, i.e. of the tape with stacked 
graphite layers, is shown in the inset of Fig. [TJ Scan- 
ning across the sample in the micro-Raman experiment, 
see Fig. [Tl we only detected the signal typical of bulk 
graphite 19|. This signal corresponds to the light areas 
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FIG. 2: (color online) Transmission spectra at selected mag- 
netic fields in the low (a),(c) and high (b) energy windows. 
The lines denoted by individual letters exhibit a \/B-scaled 
blue shift and are related to the H point. The low-energy fea- 
tures in (c), nearly linear in B, originate from the K point. 
The black arrows in (c) denote the most developed transmis- 
sion minima whose positions were plotted into the inset of 
Fig. E] For clarity, successive spectra in (a),(b) and (c) are 
shifted vertically by 0.05, 0.15 and 0.15, respectively. 

in the picture, whereas the other dark areas (tape) ex- 
hibited no graphite or graphene signal. The layers of the 
bulk graphite of various thickness cover roughly 50% of 
the tape surface. The transmission experiment was per- 
formed on a macroscopic round-shaped sample having 
5 mm in diameter. The average thickness of the graphite 
layers, ~ 100 nm, was roughly estimated from the trans- 
missivity of the sample in the visible range. 

The FIR experiments have been performed using the 
experimental setup described in Ref. [20| . To measure 
the transmittance of the sample, the radiation of globar, 
delivered via light-pipe optics to the sample and detected 
by a Si bolometer, placed directly below the sample and 
cooled down to a temperature of 2 K, was analyzed by a 
Fourier transform spectrometer. All measurements were 
performed in the Faraday configuration with the mag- 
netic field applied along the c-axis of the sample. All the 
spectra were taken with non-polarized light in the spec- 
tral range of 10-300 meV, limited further by several re- 
gions of low tape transmissivity (see gray areas in Figs. [2] 
and [3]). The transmission spectra were normalized by the 
transmission of the tape and by the zero-field transmis- 
sion, thus correcting for magnetic field induced variations 
in the response of the bolometer. 

Typical spectra, characteristic of a number of samples 
prepared in the same way, are depicted in Fig. O Sev- 
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FIG. 3: (color online) Positions of absorption lines related to 
the H point as a function of \/B. The dashed lines represent a 
least squares fit. Inset: K point related transmission minima 
showing a nearly linear with B dependence, as compared to 
the most prominent transitions (solid lines with slopes 2.3, 
4.1, 7.3 and 9.1 meV/T) anticipated after Ref. [i]. The lowest 
observed minimum is broad and can be composed of two lines. 

eral absorption lines showing a -N/B-dependcnce of their 
energies are observed in the spectra, as seen in Fig. [S] 
We relate these transitions to the H point of graphite. 
The dominant lines in the spectra are denoted by capital 
letters. The energies of the subsequent lines B to E scale 
as 1 : + : + : (2-h\/3). The integral inten- 
sities of these lines increase with the magnetic field. The 
widths of the absorption lines also increase, see e.g. the 
C line in Fig. [21 whose width increases from « 10 meV at 
B = 1.3 T to more than 30 mcV at B = 9.0 T. Weaker 
intensity features, which form a second series of absorp- 
tion lines, are denoted by Greek letters. In addition, 
features which shift linearly with _B, originating from the 
K point of graphite, are also observed in our spectra, 
see Fig. [51:;. These features, which have already been re- 
ported in magneto-reflectivity data are not further 
discussed in this work. The weak intensity modulation 
of spectra (visible in Fig. [2^,c above B k, 0.5 T), which 
does not shift with changing magnetic field, is attributed 
to interference effects in the sample. 

To interpret \/^-dependent features in our FIR spec- 
tra, we sketch the simplified SWM model of Landau levels 
(LLs) in the vicinity of the H point [^, [l3| . Starting with 
four TT-bands, Ei,E2 and the doubly degenerate E^, we 
obtain four LLs for each index n > 1 at a finite mag- 
netic field B. In addition, we obtain three levels and one 
level for n = and n = — 1, respectively. In the follow- 
ing, we use the notation i?", EJ^ for n > 0, E^^,E^_ for 
n> 1, E^ and E^^ to emphasize the band profile of the 
individual LLs. To distinguish among levels in Fig. [31 
gray color is used for LLs i?" and £'2 . The assumption 
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of zero trigonal warping (73 = 0) @, simplifies the 
eigenvalue problem to the diagonalization of 4 x 4 ma- 
trices. At = 0.5, we obtain the analytical solution 
E^^ = A, E-^ = and for n>l: 
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(1) 



where "+" and "-" correspond to E'^^'^ and E'^^'^, 
respectively, and ^ is related to 70 via the expression 
^ = 3j^eay{2h), with oo = 0.246 nm [2l|. 

The selection rules for dipole-allowed interband tran- 
sitions at kz — 0.5 predict the absorption lines at ener- 
gies (n > 1): 
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which correspond to transitions i?3_ 
with E'^-^- 



, 77111+1 



together 

and, due to the expected electron- 
hole symmetry, also to transitions E^^^^E^, and 



E^- 



3+ 

In addition to interband transitions, two 



dipole-allowed intraband transitions E'^^-^s ^ (^"^-i) 
and E^_-^E^ (^ijJo), split in energy by A, 




-3- E 

A = -5.0 meV 
y =3150 meV 



hLOo^hio-i + A=j + ^^+^B, (3) 

were experimentally found and discussed in Ref. [icH . 

A number of dipole-allowed transitions expected at 
kz ~ 0.5 are shown on the right-hand side of Fig.|4]by ver- 
tical solid arrows. Following the general assumption that 
the Fermi energy at the H point is negative, in Fig. [5] we 
assign the calculated lines to the experimentally observed 
lines, using the same notation by capital letters B,C,D,E 
and F as in Fig. [2l It is important to note that a non- 
zero value for the parameter A implies the splitting of the 
strongest B transition into two Bq and B_i components. 
However, this splitting is not resolved in the experiment. 
The upper limit of the pseudogap, |A| < 10 meV, es- 
timated from the measured width of the B line is in a 
good agreement with Ref. but in an apparent con- 
tradiction to Ref. and recent theoretical calculations 
in Ref. In Fig.gl we have used A = -5.0 meV. The 
clear observation of the B line, which accurately scales 
as \/B down to B w 0.3 T, allows us to estimate a Fermi 
energy, Ep « —20 meV, in accordance with generally 
accepted values, see e.g [l3|, but somewhat smaller than 
the ARPES result {Ep « -50 meV) 

Another important feature of our data is the clearly 
asymmetric shape of the B line as well as the significant 
line broadening with increasing magnetic field. We pre- 
sume that both these effects indicate some contribution of 
the transitions at kz ^ 0.5 to the observed transmission 
spectra. At low magnetic field, the fc^-dispersion in the 
vicinity of the H point is steeper for E2 than for E^_ and established 
this qualitatively explains the low field asymmetry. The 
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FIG. 4: (color online) LLs in graphite from n = —1 to 7 as a 
function of momentum kz in the vicinity of the H point {kz = 
0.5) at _B = 1.0 T calculated using the SWM model neglecting 
trigonal warping (73 — 0) . The arrows on the right-hand side 
show transitions considered in our interpretation. Except for 
70 and A, the TB parameters used in the calculation were 
taken from Ref. [13i |. 



dispersions become more symmetric with increasing B, as 
does the line shape of the B line. Nevertheless, the dis- 
persions of all levels are steeper at higher magnetic fields, 
which qualitatively accounts for the broadening of the 
absorption lines. However, the agreement remains only 
qualitative. An additional mechanism for the line broad- 
ening is electron scattering probability, which increases 
linearly when moving away from the Dirac point [23 |. 
Summing over the transitions for different momenta kz, 
in terms of a simple calculation of the joint density of 
states, does not help to reproduce the spectral shapes of 
the observed lines. Some additional kz dependent selec- 
tion rules, enhancing the oscillator strength in the close 
vicinity of kz = 0.5 point need to be introduced. 

The relatively small value of the pseudogap A is con- 
sistent with the very good linearity of data with -s/B (see 
Fig. [3|) . The limit of A — > allows us to rewrite Eq. ([1]) 
as E'g^ — E^2^ — ±y/^Bn, which has the form of LLs 



in a 2D system of massless Dirac fermions. The energies 
of the dipole-allowed transitions ([2]), ([3]) can then be for 
n > written as huj„ 



— ^/^{y/n + Vn+T) and a clear 
correspondence to the optical transitions in graphene is 
2^, 2^. The Fermi velocity can then be ex- 
pressed in the form c = y/^J{2eK) and from our data 
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evaluated to be c = (1.02 ± 0.02) x 10^ m/s, which is 
about 10% higher than in lisll but in a good agreement 
with the most recent data |l6l . [l7j . We stress that the 
entire fan chart of the observed inter LL transitions re- 
ported here can be described with a single parameter 
(Fermi velocity). This is a surprising observation, as 
Kohn's theorem is not expected to hold in a system with 
a strongly non parabolic one particle dispersion, so that 
electron-electron interaction may differently alter excita- 
tions between different pairs of LLs. Note that the linear 
fits in Fig. [3] extrapolate to an onset ~ 3 meV instead of 
zero, which possibly suggests, apart from the anticross- 
ing splitting A, an additional mutual shift of the electron 
and hole Dirac cones, however, this effect lies on the edge 
of the experimental accuracy. 

We turn now our attention to the additional transi- 
tions, denoted in Figs. [2] and [3] by Greek letters, which 
are in general of a lower intensity. They exhibit a clear 
-v/B-dependence and have no counterpart in the trans- 
mission spectra of graphene [2^ . To interpret these 
lines we have to go beyond the simplified model used in 
0, [13] and consider a more complete analysis [2^ . The 
selection rules n— >n ± 1 (rt > 1) is predicted not only 



for the main transitions i?5 



^E^^^ and 



£;}'±\ but 



also for a second series of lines, presumably weaker in in- 



tensity, E^_ 



^E'l^^ and EJ^- 



^E. 



Tl±l 

3+ 



We can then relate 



the absorption lines a,j,S and e to transitions which 
are symmetric with respect to the Dirac point, as de- 
picted in Fig. [4] by the vertical dotted arrows. In agree- 
ment with expectations, the energies of these lines scale 
as 1 : : : 2 and the same Fermi velocity is derived 
(c = (1.02 ±0.02) X 10^ m/s) as for the main transitions. 

Another series of weak absorption lines, predicted 
in [2^, should satisfy the selection rules n^n± 2. Their 
oscillator strength is directly connected to the parame- 
ter 73. The /3 line can then be assigned to transitions 



e: 



3(1) 



iSgiJ*'' and/or E^^'^' ^E'^'^^' , representing thus a 
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HOPG can, with a reasonable accuracy, be described as a 
purely 2D gas of Dirac fermions. This justifies the model 
adopted for interpreting STS experiments [13] and quan- 
tum oscillations in graphite [ll|, Nevertheless, the 
latter analysis is not consistent with our data, when es- 
timating the hole density. Whereas in our experiment, 
the hole filling factor v — Q is achieved at _B « 0.3 T, 
when the B line clearly appears in the spectrum, the 
same value of v is found at ma gnet ic fields above 1 T in 
the magnetotransport data ll|, |l^ . 

On the other hand, the series of additional absorp- 
tion lines of weaker intensity, which are dipole-forbidden 
in the strictly 2D gas of Dirac particles, show that in 



HOPG we are dealing with a strongly anisotropic but 
nevertheless 3D system and the model of a purely 2D 
gas is not valid in general. The 3D nature of graphite is 
also shown by the pronounced low field asymmetry of the 
B line and possibly by the observed increase in linewidth 
with increasing B. 

In summary, FIR magneto-transmission measurements 
have been used to probe the Dirac holes at the H point 
in a thin graphite layer prepared by the exfoliation of 
HOPG. We find a relatively small value for the pseudo- 
gap ]A] < 10 meV, consistent with the observation of 
Dirac fermions at the H point of graphite, with a Fermi 
velocity of c = (1.02 ±0.02) x 10^ m/s. The main absorp- 
tion lines can be understood using a model which assumes 
a strictly 2D gas of Dirac fermions, giving evidence for 
the strong anisotropy of HOPG. Nevertheless, the pres- 
ence of additional weaker transitions, dipole-forbidden in 
a purely 2D system of Dirac fermions, can only be un- 
derstood by considering the 3D nature of graphite. 

The present work was supported by the Euro- 
pean Commission through Grant No. RITA-CT-2003- 
505474, by contract ANR-06-NANO-019 and projects 
MSM0021620834 and KAN400100652. 



direct indication of the trigonal warping in graphite, 
demonstrated experimentally e.g. in [ISj, [l6|. 

Having interpreted all observed absorption lines, we jg 
can draw the following conclusions. The dominant tran- [7" 
sitions in the spectra, which have their counterpart in 
the spectra of graphene [l^ 24], justify that holes in [8 
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